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Abstract 

We introduce a simple parametrization of the dark energy equation of state, u;, which is 
motivated from theory and recent experimental data. The theory is related to the tracker 
solution to alleviate the fine tuning problem. Recent experimental data indicates that the 
present value of uj is close to —1. We analyze the evolution of uj from the separation 
of CMB peaks and the time variation of the fine structure constant, a. We find that 
-1.00 < cu(°) < -0.9711^27 and 1.76lg:i x 10"^ < {duj/dz),=o < QMltl^^ at 95% 
confidence level (CL). 



1 Introduction 



Recent cosmological observations show that a component with negative pressure (dark 
energy) should be added to the matter component to give the critical density today [1]. 
The cosmological constant and/or a quintessence field are the most commonly accepted 
candidates for dark energy. The latter is a dynamical scalar field leading to a time dependent 
equation of state, uj. 

There are certain hopes that quintessencelike models may help alleviate the severe fine- 
tuning problem associated with the cosmological constant and that the nature of dark 
energy can be understood by measuring the behavior of uj with respect to time [2, 3]. But 
it may not be practical to test every single quintessence model by using experimental data. 
So a model independent approach for quintessence could be an effective way to study the 
properties of dark energy. 

To investigate the cosmological evolution of a;, we may use the relation between the 
luminosity distance and redshift {z) , which can be obtained from type la supernovae (SNe 
la) as standard candles [4]. In this case, the distribution of the supernovae (SNe) along z 
provides a clue for the determination of uj. The dependency of the luminosity distance on oo 
increases in proportion to z. The separation of cosmic microwave background (CMB) peaks 
can also be used to investigate uj because the location of peaks depends on the amount of 
dark energy today and at last scattering as well as u [5]. In addition, the cosmological 
evolution of the fine structure constant, a may be used to determine the z dependence of cu 
[6] . The advantage of using a is that there will not be any degeneracies between parameters 
due to the multiple integrals as in the determination of luminosity distance [7] . 

A number of parametrizations of cu have been used to probe the nature of dark energy. 
The simplest is the averaged value Cu, defined as cu — f uj{z)ri^{z)dz/ J D,^{z)dz where fi^ 
is dark energy density contrast. This parametrization is weighted by indicating that 
u is more significant if is a larger part of the critical density. If uj does not change 
during the recent history of the universe, when is comparable to matter density (l^m), 
then the average uj'^'^^ will be equal to present value, uj'^'^\ The next intuitive method is a 
first order Taylor expansion of uj in z, uj{z) = uj\z=o + du/dz\z=oz. This diverges at large 
value of z, which makes it only useful for low z analyses, such as SNe [8]. Another simple 
parametrization is a; = ujo+ujiz/{1 + zY where p — 1,2. When p is equal to one, ujq+uji is 
constrained as a negative value to satisfy SNe data. In addition, at decoupling should 
be neghgible based on WMAP [9]. However with this constraint we can not get a small 
value of Q(f, at decoupling. The p = 2 case can solve this problem [10]. A simple extension 
of the previous cases is given by uj{a) = uJoUJi{a'^ + a^) / {ujia'^ + ujoaj), which includes huge 
fiexibility and the intuitive role of each parameter [8]. In this case, both ujq and uji should 
have the same (negative) signs. Otherwise, uj will diverge at = — (a;o/to'i)a^. There have 
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also been more complicated parametrizations [11]. 

We introduce a new parametrization of cu, which looks similar to that in [8]. However 
it is necessary that the sign of cuq is opposite to that of cui to have a smoothly changing cu. 
We will not consider the case of a; < —1, which can be obtained in scalar-tensor gravity 
models [12], phantom models [13], and brane models [14]. 

This paper is organized as follows. In the next section we introduce a parametrization of 
ou and specify some of its parameters. By changing these parameters we can mimic several 
quintessence models. In section 3, we use CMB peaks to restrict the parameters. We check 
the time variation of a based on this parametrization in section 4. Our conclusion is in the 
last section. 



2 Parametrization of uu 

For cosmological evolution equations, it is convenient to introduce a variable x as the 
logarithm of the scale factor a, 

x = lna = -ln(l + z) (2.1) 

where we choose the present scale factor a^^^ = 1. We propose a simple parametrization of 
the dark energy equation of state, 

Uj{a)^UJr 5"^ =^ uixj^UJr 7 ^ 7 ^ (2.2) 

ai + at exp(gx) + e'xjgyqxc) 

where (equally Xc) , and q are constants. Instead of leaving uir and uiq as arbitrary 

constants, we adopt the tracking condition in the early universe i.e. cu changes as that 
of radiation (uj = ujj. = 1/3). Also we use the experimental evidence which indicates cu 
approaches negative one at present. These constraints fix two constants in Eq. (2.2) : 

1 , , 

^r- =3, t^o = -3 (2.3) 

There remain two arbitrary constants Oc and q, which describe the scale factor at 
changeover and the duration of it respectively. When a is close to Oc, ui starts to de- 
crease and cu reaches to —1 as a approaches a° = 1. We call a transition of uj from the 
tracking region to the c<j ~ — 1 region as changeover. We can also narrow the range of the 
Qc value, when we consider the fact that initial value of is expected to be 120 orders of 
magnitude larger than the present value. We will see this in the next section. As we can 
see in Eq. (2.2), for larger value of q, the duration of changeover of cu decreases. During 
early times a; ~ a;^ and it stays with same value before a approaches Cc. cu will approach 
cUr X cuq when a ^ Oc. This is the reason we choose the values in (2.3). If Oc is too close to 
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Figure 1: The evolution of uo. (a) The dependence on the changeover scale, Xc of the 
evolution of a;. A smaller value of Xc means an earlier change of uj shown as the dashed 
line, (b) The dependence on q of the evolution of uo. The solid line (g = 1) shows a slower 
change than the dashed one (g = 4). That means the duration of the changeover becomes 
longer as the value of q becomes smaller. 

a*^°) = 1, then the present value of uj will not be close to —1 for small q. Thus there is a 
limit to Xc for each value of q. The cosmological evolution of uo due to this parametrization 
for different values of q and Xc are shown in Fig. 1 and Fig. 2. The effect of Xc on the 
evolution of uj is shown in Fig. 1 (a). A smaller value of Xc gives an earlier change of uj. As 
Xc decreases, the period of = — 1 increases. A smaller value of q shows a slower change of 
UJ as in Fig. 1(b). For a specific value of Xc = —2.64, we require q > 0.62 to satisfy WMAP 
data. For larger values of q we can still get uj^^^ = — 1 even with large Xc (i.e. the duration 
of changeover is small). From Fig. 2 we show the minimum value of q [qmm) for each Xc to 
satisfy WMAP data (u;*-'^'' < —0.78). For a smaller value of Xc we allow a smaller value of 
Qmin to satisfy WMAP due to the fact that u changes early and it takes more time for uj 
to reach today when Xc is small. This is shown in Fig. 2(a). We also show the restriction 
of q for a specific value of Xc in Fig. 2(b). The advantages of this parametrization are its 
simplicity and the possibility of mimicking various models. However this parametrization 
cannot account for an oscillating uj, which can be obtained from a cosh(0)-potential and 
the coupled cases with general potentials [15, 16]. Thus this parametrization is good only 
for monotonically varying uj. 
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Figure 2: The present value of the equation of state, u^'^\ (a) gmin means the minimum 
value of q for each different value of Xc to satisfy WMAP data (0;^°-' < —0.78) [9]. (b) The 
shaded region is forbidden from WMAP. From this we can make the restriction q > 0.62 
when Xc = —2.64. The reason for the choice of this specific value of Xc will be shown later. 

3 CMB peak spacing 

The location of the acoustic peaks and the spacing between the peaks can be estimated if 
an adiabatic initial condition and a flat universe are assumed [5, 17]. The acoustic scale, 
Ia is given by the simple formula 

La = TT — = vr = vr— (3.1) 

S S CsTis 

where (Ia is the angular size distance to decoupling, tq and tis are the conformal time today 
and at last scattering. The sound horizon at decoupling is given by s = CsTis and the 
average sound speed before last scattering is, 

c^ = It^ where c- = 3 + ?4t (3.2) 

where pb and pr are the baryon and photon energy densities, respectively. The location of 
the m-th peak and trough is slightly shifted by driving effects and this can be written as 

[18]; 

Im = Uim - ipra) = - Ifi - 5^rn) (3.3) 

where = y^i is the overall peak shift and 5iprn is the shift of the m-th peak relative to the 
first. We attach the fitting formulae in the appendix for the sake of completeness. To see 
the effect of quintessence, we start from the Friedmann equation. 

where if is the Hubble expansion rate, = M^/S TT is a reduced Planck mass, and pi are 
the energy densities of each component. From the continuity equation of energy density we 
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find, 

d\np^ — —3(1 + u!)dx. 
This equation can be integrated using the parametrization of cu in Eq. (2.2) : 



(0) 



i(a^) = P<t> 



l + exp(,..) j 



(3.5) 



(3.6) 



where p^^^ is the dark energy density today. We also find the dark energy density contrast 
(f2<^) from the above equation (3.6). 



Per 



1 + -^{a + tteq) 



(0) 



1 + aI 



(3.7) 



where Q-'^-* = pf^ / p^^r present energy density contrast of each component. In the 

second equality we use the relation Vtf^ = VL^agq where agq is the scale factor when the 
radiation and matter densities are equal. 



Table 1: Cosmological parameters used in the analysis. We use WMAP data [9]. 









Zls 






0.73 


0.27 


0.0224 


1089 


3233 


0.93 



Prom Eq. (3.6), we can see that the dark energy density, evolves as the energy density 
of radiation {pr oc a~^) during a <S ac- Ptj, starts to departure from this behavior as a 
approaches Qc- This explains the tracking behavior of dark energy during the radiation 
dominated era. The cosmological evolution of the dark energy density is shown in Fig. 3. 
p^ is changed by 120 orders of magnitude through the entire history of the universe. Also 
from Eq. (3.6) we know that the initial value of dark energy, p'^^\ac/ cnY is independent 
of g, where aj is the scale factor at the beginning of the universe. If we use the fact 
that p\l p^^ ~ 10^^° and Oj ~ 10~^^, then we find that Oc should be approximately 10~^ 
(equally Xc ~ —4). We can make further restriction on Oc from the following consideration. 
To be compatible with observational data, the energy density of quintessence must be 
subdominant during Big Bang Nucleosynthesis (BBN) [3], Q^^^^\x ~ —23) < 0.2 at 
T ~ 1 MeV. If wc consider the fact that 1 3> Oc ^ ^ clbbn-i then the equation of dark 
energy density contrast (3.7) is approximately : 
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Figure 3: (a) The cosmo logical evolution of for specific q and Xc values. For other 
values of these parameters, we can get similar behaviors which can account for 120 orders 
of magnitude change in p^. (b) The evolution of when we fix = 0.73. When 
Xc = —2.64, a smaller value of g = 1 (solid line) shows a slower change in fi,^ at late 
universe compared to g = 4 (dashed line). For a smaller value of Xc = —4.00 (dotted line), 
fl^ is negligible during early universe for any value of q. 



^(t>{cLBBN) 



1 _|_ "1 "eg 



< 0.2 ^ a, ^ ( ^(o/ ^_'^(BBiv) < 0.073 (3.8) 



This corresponds to Xc < —2.60 independent of q. As previously mentioned, the slope of 
during the radiation dominated era has nothing to do with Oc or q. The effect of q appears 
when a approaches Oc. This effect on p^ is hardly seen as in Fig. 3(a) due to degeneracy. 
However we can see this effect from Fig. 3(b). A larger value of q (dashed line) gives the 
steeper change in p^ and than those of a smaller value of q (solid line) at late time. As 
we can see in the equation (3.8), the Q^j, value depends only on Oc before a reaches to aeq. 
Thus for the smaller value of a^, value will be more close to the value of cosmological 
constant, Qa. The evolution of dark energy density shows degeneracy between the different 
choice of parameters on uj. This may require additional experimental data to investigate 
the evolution of u instead of using SNe. Because SNe data is analyzed by the luminosity 
distance, which mainly depends on p^. In the above graphs we fix fi^^-* = 0.73. Now we 
rewrite the equation (3.4) by using conformal time, r : 



(l-0<»)-f!f)(. + a.,)+Of(^Li|)'"| (3.9) 

where Hq is the present value of Hubble parameter. If we use WMAP data ~ 4.9x 10-^ 
fiW ^ 0.27, n'f^ ~ 0.73, and a^g = 1/(1 + z^q) = 1/3234, then we have the following 
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X =-2.64 

c 



1.5 2.0 2.5 3.0 3.5 4.0 

q 

Figure 4: The locations of the first and the second CMB peaks. The uptriangles indicate 
the first peaks, h and the downtriangles show the second peaks, I2 as a function of q. Both 
h and I2 converge as q increases. 



approximation of the above equation (3.9). 




(3.10) 

The approximation comes from the fact that Q^^^ ^ ^m,4>- indicate the inequahty 
between the first and the second expressions, which is frequently used in the literature. 
The numerical difference between two expressions is not negligible for the specific values of 
Gc and q. So we will use the first (accurate) expression in the following calculations. From 
this equation we can find the numerical value of r when ac and q are fixed. 

dr = H,^ ((1 - <))(a + ae,) + {"^f) ^^'^ (3.11) 



Table 2: The spacing and the location of the CMB peaks for several values of q when 
Xc = —2.64. Here we use the cosmological parameters in Table 1. 





q 


Ia 


h 


k 




-2.64 


1.5 


300.8 


218.3 


535.0 


5.17 




2 


302.4 


219.3 


537.7 


5.36 




3 


303.1 


219.8 


538.9 


5.40 




4 


303.1 


219.8 


539.0 


5.41 



By using the above equations we can find the acoustic scale {I a), the locations of the 
first two peaks {hyh), and at last scattering {^''I). They are given in Table 2 and Fig. 
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Figure 5: The 68 % (light shaded) and 95 % (medium shaded) confidence allowed regions for 
Xc and q using WMAP data. Asterisk indicates the best fit values [xc = —2.64, q = 4.01). 
The dark shaded regions are out of 2a confidence level. The constraint from the dark energy 
density at BBN determines the upper limit on Xc (< —2.60). The lower limit of Xc (> —2.72) 
is obtained from x^-analysis. There is no upper limit on q values from WMAP. However 
we can obtain the upper limit on q from the analysis of the time varying a in the following 
section. 

4. For a value of Xc, the location of CMB peaks depends on g-values. However CMB peak 
values converge to some values as q increases. From WMAP measurements the acoustic 
scale and the locations of the first two peaks are [9] : 



If we assume that the errors are Gaussian and uncorrelated, then we can find a x^-statistic 



where li^obs is the observational data of each peak, h^theory is the theoretical value of /j, and 
(Tj is the statistical uncertainty for each peak. The likelihoods for the parameters (g, Xc) 
of the parametrization given in (2.2) are calculated from so that the 68% and 95% 
confidence regions are determined by Ax^ = ~ Xo = 2.30 and 6.17 respectively. Xo is 
for the best fit model found. The best fit values to the WMAP data is obtained from 
Xc = —2.64 and q = 4.01. The likelihood contours for WMAP are shown in Fig. 5. The 
68 % (light shaded) and 95 % (medium shaded) confidence allowed regions for Xc and q are 
shown using WMAP data. There is no upper limit on q values for several Xc values. From 
this we find at 95 % CL : 



/a = 301±1, /i = 220.1 ± 0.8, /2 = 546 ±10. 



(3.12) 




(3.13) 




(3.14) 
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cu is almost a constant at present, which is consistent with recent observation [19]. 



4 Time variation of the fine structure constant 

In this section we consider the 0-dependence of the gauge couphngs in Standard Model 
(SM), Bf{4>) [16]. The interaction of a light scalar field with electromagnetic field is 
represented in the Lagrangian as Cem — — -Bf(0)-F^;/-F'^'^/4, where F^^^ is the electromag- 
netic field-strength tensor. The main assumption in this work is that all functions Bf{4>) 
and V{(1)) admit a common cxtremum, which is a generalization of Damour-Nordtvedt and 
Damour-Polyakov constructions [20] . To this end, we proposed an ansatz that would relate 
BF{(t)) and V{(i)). 

Bf{<P) = ( ^-£±Xm^y\ where 6^ + 1 > (4.1) 
V 1 + Of / 

where np and bp are two dimensionless parameters allowing us to cover a wide range of 
possibilities. We will assume that dark energy is a quintessence field and its potential is 
parameterized by cu. So instead of introducing a potential of the scalar field, we will use 
cu. We can rewrite (4.1) by using the parametrization of uj (2.2). We also use the relation 
between the potential V{(f)) and the energy density of the scalar field, 

Vi'P) = ^ V{x) = pf exp(-4a:)^i(x)^2(x) (4.2) 



where 



( exp(gx) + cxp(gxe) y^^ , , 

= [ l+exp(,.J j ■ <^-^' 

^2^^-) ^ 3exp(gx) +exp(g3;,) ^^^^ 
3^exp(ga;) -|- exp(ga;c'^ 



With these equations (2.2), (3.6), and (4.2), we can rewrite the gauge coupling in Eq (4.1) 
as follow, 

B,(.) = ('V + ''x„(-LO^.M.4.(..)/.4.(0)j'- (^^^j 

where we assume that Vq — V{x — 0) to satisfy Bf{Qi) — 1. We can find the time variation 
of the fine structure constant by using the following relation between Bp and a, 

Aa ^ a(a;) -a(0) ^ 1 _ ^ 

a q;(0) Bf{x) ' ^ ' ' 
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Figure 6: A typical evolution Aa/a for specific np, bp, q and Xc values. We can get similar 
behaviors with other values of parameters, np and bp values are obtained from the result 
of Murphy et al. Several values of Act / a for various choices of the parameters at different 
z are indicated in Table 3. 

From equations (4.5) and (4.6), we find 

Aa _ / l + bp 

a ~ + exp(-4x)Ai(x)A2(x)/A2(0) 




Table 3: The values of Aa/a at the different z for the several values of q when Xc = —2.64. 
Here we use the QSOs result to fix other values [21]. np has been scaled by a factor of 10^, 
Aa/a values have been scaled by a factor of 10^ except the value at decoupling {Aa/a)is, 
which is given in the last column without any scaling. 





Q 


Hp 


bp 


( A.a\ 


I — Jl.5 


( Aoi\ 




I — j 1089 


-2.64 


3 


4.25 


11 


-5.4 


-1.25 


-0.017 


-0.041 


-0.057 






0.71 


1 


-5.4 


-1.26 


-0.017 


-0.042 


-0.011 




3.5 


3.68 


3 


-5.4 


-1.01 


-0.023 


-0.025 


-0.053 






1.84 


1 


-5.4 


-1.01 


-0.023 


-0.025 


-0.028 




3.9 


4.01 


1 


-5.4 


-0.84 


-0.016 


-0.016 


-0.060 



In Table 3, we display Aa/a values at different redshifts for the various parameters. Es- 
pecially we show specific results for Xc = —2.64, because this gives the best-fit value for 
the separation of CMB peaks as explained in the previous section. Aa/a at the cosmic 
microwave background epoch, which is constrained as {Aa/a)is > —0.06 [22] are shown in 
the last column. There are four parameters {np,bp,Xc, and q) to be fixed. In Fig. 6 we 
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Figure 7: The cosmological evolution of a for the various parameters. Murphy et al data 
(medium shaded) is used in all of the figures, {^a/a)z=z = (—0.543 ± 0.116) x 10~^. We 
use Xc = —2.64 for all of the graphs except (d). (a) The evolution of Aa/a with respect 
to bp > 1 for Up = 5 X 10^^. When g = 3 (dotted line) a changes slower than when q = 4 
(solid line). Thus we have more acceptable values of b for the smaller q (10.6 < bp < 16.9 
when q = 3). (b) The effect of changing np on Aa/a while we fix bp = 10. (c) and (d) the 
effect of varying q and Xc on the variation of a, respectively. 
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show a typical evolution of Act/o; for a specific set of four parameters {i.e. np — 4: x 10"^, 
bp — 1, q — 3.9, and Xc — —2.64). We can change four parameters to explain the effect 
on the evolution of a. However it is impossible to fix all four parameters at once. Thus 
we show the effect of the changing each parameter with fixing all other parameters at each 
time. We sec that smaller values of q lead to more room for the parameter space regions 
of np and bp, because smaller values of q lead to slower change in a. For given g-values, 
there may not be any solution to fit experimental data as shown in the first panel of Fig. 
7. However as we increase the feir-values we can satisfy the experimental data. We have 
larger bp range for smaller values of q. Similar consideration can be made for varying np. 
For example, even though the dotted lines (q = 3) in the second panel does not fit experi- 
mental data when np > 3 x 10~^, as we decrease np we can recover the suitable Aa/a for 
given q and Xc. The medium shaded regions of each panel are allowed from Murphy et al, 
{Aa/a)z=3 = (— 0.543±0.116) x 10~^ [21]. The effect of changing q and Xc are shown in (c) 
and (d) of Fig. 7. As we expect there is lots of room for the choice of q and Xc because we 
have extra degrees of freedom {up and bp) to be fixed in addition to the original parameters 
of uj. However combined with the analysis of the separation of CMB peaks, the parameters 
on gauge coupling are much narrowed as in Table 3. We can make the allowed regions for 
np and bp from the narrowed parameter spaces of Oc and q. This is indicated in Fig. 8. The 
sloped lines are obtained from the normalization of Murphy et al data for the given q values 
when Xc is fixed as —2.64. The vertical lines for each q come from both the CMB bound 
on a and the assumption of fo > 1. There is only one allowed value bp = 1 for q = 3.9 ( 
solid line ). This gives the upper limit on g < —3.9, which we can use for the consideration 
of CMB peaks. For q = 3 (dotted line), there are more parameter spaces to satisfy the 
constraints. In this case bp can be increased up to 11. Obviously this gives much narrower 
parameter spaces for each parameter compared to Fig. 7. Now we can restrict {np, bp) ~ 
(10-^ 1) as in [16] and 1.5 < q < 3.9 for Xc = -2.64. 



5 Conclusions 

We have introduced a simple parametrization of cu based on both theory and experiment. 
A tracking solution is proposed to alleviate the fine tuning problem and the present value 
of cu is observed to be close to —1. We have analyzed the separation of CMB peaks and the 
time variation of the fine structure constant with this parametrization. We have obtained 
the constraints of the parameters of uj from these analyses which use both low-redshifts and 
high-redshifts data. We have been able to use the analytically integrated value of the dark 
energy density contrast, from this parametrization. 
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bp 

Figure 8: The constraints on np and bp from both the result of Murphy et al [21] and that 
of CMB [22]. np values in the graph are scaled by 10^. The sloped lines are obtained from 
QSOs result. The vertical line for each case is obtained from CMB. The shaded regions 
are allowed ones from both constraints. There is only one 6i?-value for q = 3.9 (solid line) 
when we assume that 6 > 1. So we can get the upper limit g < 4. 



From WMAP data, we have obtained a very specific value of Xc ~ —2.64 (equally ~ 
13) for the best fit value. When Xc = —2.64 we have obtained the very small {dw / dz) ^=0, 
which is consistent with the recent observations, (dw/dz) depends on Xc, which is related 
to the epoch when the time evolution of u begins. A smaller value of Xc causes an earlier 
evolution of u compared to a larger value of Xc- With a smaller Xc, evolves to its present 
value p^^'' at an early epoch and it stays with this value because ci; ~ — 1. In this case 
VL^ is negligible until it reaches its present value and that is similar to when we have the 
cosmological constant as dark energy. However for some values of Xc the value of Vt^ can 
be significant at early times. This behavior is quite different from that of a cosmological 
constant. Thus the difference between the cosmological constant and quintessence is obvious 
for some values of Xc- Due to CMB consideration we have focused on x^ = —2.64 and this 
has shown the big discrepancy between quintessence and the cosmological constant. We 
have also found that to be consistent with BBN {i.e. Vt^^^^"^ < 0.18) we need q > 1.50 
when Xc = —2.64. Thus this may give one clue for the nature of dark energy. We have 
made the restriction to the value of q, q^i^ from the constraint on Xc- As Xc increases, we 
will get bigger gmin as in Fig. 2. We have also found the upper limit on Xc when we consider 

( Ft Ft /VI 

the limit of , Xc < —2.60. We have reached to the lower limit q > 1.50 when we used 

~ 0.01 from WMAP. We have had the similar lower limit on q from the x^-statistic. All 
of this consideration gives the unique lower limit of g > 1.50 when we choose Xc = —2.64 
However there was no upper limit of q from the CMB analysis alone. 

We have introduced two more parameters [np, bp) from the generic form of gauge cou- 
pling to analyze the evolution of Aa/a. However we had only three parameters to be fixed 
after specifying Xc- We have found the viable parameter regions for each parameter by 
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using Ack/q; data. Prom this we have found the well known parameter regions of gauge 
coupling {nF^hp) ~ (10~^, 1) [16]. We have also found the upper limit q < 3.90 with this 
analysis. 

We have found a viable parametrization of uj with the very specific parameter regions, 
Xc = —2.64 and 1.50 < q < 3.90 from the analyses of both the separation of CMB peaks 
and the time variation of a. We have also found that —1.00 < cu^^^ < — 0.971]'^o;o27 and 
1.76l[J:i X 10-^ < {dLj/dz),=o < OMltomi at 95% CL. Thus we may not see any change 
of u at present. Even though the present values of uj and duo/dz are almost identical to the 
values of the cosmological constant, VL^ evolves from BBN to the present with significant 
values. Thus this model has the different cosmological behaviors compared to ACDM 
model. We may need to check the similar parametrization and/or more data to investigate 
the nature of dark energy. We may also need to modify uj to have more than one changeover 
which can show a non-monotonic behavior. 
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A Appendix 

We have used the analytic approximations for the phase shifts Doran, which we show here. 
The overall phase shift is given by 



(p = (1.466 - 0.466n^) + 0.291O; 



■4> 



(A.l) 



where the fitting coefficients are 



oi = 0.286 + 0.626a;b 

as = 0.1786 -6.308a;b+174.9u;^-1168u;^ 



(A.2) 



with uJi, = Q[ 'h'^, and Q''f is given by 




(A.3) 



and 



n = pr{zis) / pm{.Zls) 



(A.4) 
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is the ratio of radiation to matter at decoupling. The relative shift of the first acoustic peak 
is zero, 5(fi — 0, and the relative shifts of the second peak are given by 



S^2 = Co - cir, - C2r-^3 + 0.05(n, - 1) (A.5) 

where Ug is scalar spectral index and 

Co = -0.1 + (0.213 - 0.1230fJ X exp(-[52 - 63.61^fjwb) 

Cl = 0.015 + 0.063 exp(-3500u;^) 

C2 = 6 X 10-^ + 0.137(a;6-0.07)2 

C3 = 0.8 + 2.3Qf, + (70 - 126Qf,) u;6 (A.6) 
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